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SIGNIFICANCE  AND  EXPLANATION 

The  uniqueness  property  for  solutions  of  an  initial-value  problem  turns 
out  to  be  a  fundamental  tool  which  interacts  with  many  other  questions  of  a 
qualitative  nature*  For  an  evolution  equation  that  is  supposed  to  serve  as  a 
model  for  a  natural  phenomenon ,  it  is  often  crucial  if  one  wants  the  model  to 
be  credible r 

The  equations  of  the  type 

(E)  ufc  -  ■  0  in  (0,»)  x  , 

where  <Pi  ♦  [O,")  is  nondecreasing,  arise  in  many  applications.  These 

include  the  flow  of  a  gas  through  a  porous  medium,  the  spatial  spread  of 
biological  populations,  the  interstellar  diffusion  of  galactic  civilisations, 
the  heat  flow  in  a  material  with  a  temperature  dependent  conductivity,  the  so- 
called  Stefan-problem,  mots*  problems  in  plasma  physics,  etc.  i The  variety  of 
these  applications  explains  the  abundant  literature  concerning  uniqueness  for 
such  problems.^ 

This  paper  is  a  new  contribution  to  this  question.  It  deals  with  the 
case  when  the  initial  datum  is  a  measure,  (E)  being  understood  in  the  sense  of 
distributions.  Thanks  to  this  general  setting,  it  recovers  most  of  the 
previous  results  and  takes  into  account  a  larger  number  of  physical  models* 

But  its  main  interest  comes  from  the  fact  that  it  solves  the  *right*  unique¬ 
ness  question  in  a  form  which  arises  in  many  other  related  questions.  For 
instance  the  study  of  the  asymptotic  behavior  of  the  solutions  of  (B)  can  be 
reduced  to  this  type  of  a  uniqueness  problem  with  a  Dirac  mbs  as  initial 
datum. _ 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


UNIQUENESS  OF  THE  SOLUTIONS  OF  ufc-  A^(u)  -  0 
WITH  INITIAL  DATUM  A  MEASURE 
Michel  Pierre 

Let  us  first  state  the  main  result  of  this  paper*  Let 
<fi :  [0,*°]  +  [0 ,*°)  be  nondecreasing,  locally  Lipschitz  and  ^(0)  ■  0.  Let  us 
consider  the  problem 

|(1)ue  L^tOfT)  x  rn)  n  l"((t,t)  x  rn)  wt  e  (0,T) 

(P)  \  N 

^  (2)  ut  -  A*(u)  -  0  in  P*((0,T)  x  r") 

where  "in  P'((0,T)  x  rn>"  means  "in  the  sense  of  distributions  in 
(0,T)  x  rn». 

THEOREM  t.  Let  u  and  u  be  two  nonnegative  solutions  of  (P).  If  N  ■  1  or 
2,  assume  also  that  ^(u)  e  L^tO.T)  x  rN),  Then 

(3)  lim  ess  u(t)  -  u(t)  *0  in  P*(RN) 

t+O 

A 

implies  u  ■  u« 

A  lot  of  papers  have  already  been  concerned  with  the  uniqueness  of  the 
solutions  of  problem  (2)  especially  in  the  particular  case  of  the  porous  media 
equation  ufc  -  Au”  -  0*  See,  for  example,  Oleinik  [19],  Kalashnikov  [12], 
Gilding  and  Peletier  [10],  Kamin  [13]  for  one  space  variable,  Vo 1' pert  and 
Hudjaev  [22] ,  Sabinina  [20]  and  finally  Brezis  and  Crandall  [6]  in  the  general 
case.  In  the  latter  work  which  recovers  most  of  the  previous  uniqueness 
results  contained  in  the  above,  the  initial  value  is  assumed  to  be  in  L^R1*) 
(or  L*(*P))  and  (3)  holds  in  L1  (■?*)» 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-004 1 .  This 
material  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
Grant  No.  MCS- 7927062. 


Here  the  initial  value  is  only  a  finite  measure  which  is  the  limit  of 
u(t)  in  the  sense  of  measure  (only!)  -  its  existence  is  implied  by  (1)  and 
(2)  (see  lemma  2).  This  leads  to  a  more  sophisticated  analysis  whose  main 
difficulty  is  solved  by  using  precise  properties  of  the  potentials  of  the 
functions  u(t)  for  N  >  3.  A  different  proof  for  N  -  1,2  is  necessary  due 
to  the  non-existence  of  potentials;  it  requires  ^(u)  e  L1  which  is  in  fact 
implied  by  ( 1 ) ,  ( 2 )  in  the  cases  of  interest  ( see  remark  3  and  theorem  4 ) . 
Among  the  uniqueness  results  mentioned  above,  only  Kamin  in  [13]  considers  the 
case  of  a  measure  as  initial  data  in  the  particular  case  of  dimension  1  with  a 
Dirac  mass. 

Here  our  proof  is  quite  general  and  not  particular  to  It  would 

carry  over  to  equation  ( 2 )  in  a  bounded  domain  ft  of  H1*  with  Dirichlet  or 
Neumann  boundary  conditions. 

The  section  1  is  devoted  to  the  proof  of  theorem  1.  In  the  section  2,  we 
state  an  existence  theorem  for  solutions  of  (P)  whose  initial  value  is  a  given 
nonnegative  finite  measure.  We  also  study  the  dependence  on  the  initial  data. 

Some  last  comments  about  motivations.  The  equation  (2)  arises  in  many 
applications.  We  will  not  recall  them  here  since  they  can  be  found  in  most  of 
the  papers  mentioned  above  or  in  the  references  they  contain  (see  also  [18] 
for  a  survey  about  the  porous  media  equation) .  The  case  when  the  initial 
datum  is  a  measure  is  also  a  model  for  physical  phenomena  (see  [13]  and  [23] 
p.  677].  Moreover  it  arises  in  several  mathematical  questions.  One  example 
is  the  study  of  the  asymptotic  behavior  for  the  solutions  of  the  porous  media 
equation  which  can  be  reduced  to  the  uniqueness  problem  with  a  Dirac  mass  as 
initial  data  (see  Friedman-Kamin  [11],  Kamin  [15]). 
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SECTION  1 

We  denote  by  Cc(R^)  ( reap.  c^R**))  the  set  of  continuous  functions  on 
rf*  with  compact  support  (resp.  bounded)  and  by  #(R**)  (reap.  *J+(B **))  the 
set  of  finite  (resp.  and  nonnegative)  Radon  measures  on  R*.  A  sequence 
WncM(l^)  is  said  to  be  converging  to  it  in  o(M(RN),  Cc(R?*))  (resp. 
o(M(rN),  Cb(l^)))  if,  for  any  ye  Cc(R^)  (resp.*  e  <^(8**)) 

lim  /  *  du  -  J  „*du  . 

n—  rN  n  RN 

Case  N  >  3 

Let  E  ( x)  -  - - - — —  where  Sj.  is  the  surface  of  the  unit  N- 

N  (N-2)S  lxiN_2 

N 

sphere.  For  u,  u  solutions  of  (P)  we  denote 

a.e.  t  e  (0,T),  v(t)  -  *  u(t),  v(t)  »  ^  *  u(t)  . 

Since  u(t)  c  L^tR**)  n  L* (RN),  v(t)  e  C^R**)  n  LJ?(RIi)  for  p  >  (This 

follows  from  elementary  properties  of  the  convolution  in  it*.) 

Lemma  2.  if  u  is  a  nonnegative  solution  of  (P),  u(t)  converges  in 
c(M(R1*),  Cb(l^))  to  some  1 1  e  M+(RN)  when  t  ♦  0  essentially.  Moreover 
when  t  decreases  to  0,  v(t,x)  increases  to  v(0,x)  -  (E^  *  U)(x)  for 
all  x  c  it*. 

(Note  that  v(0,x)  is  lower  semi- continuous  and  is  the  usual  potential 
of  the  finite  measure  u ) . 

Proof  of  lemma  2 

The  relation  (2)  implies  that 


3- 


a.«.  0  <  s  4  t  4  T,  u( t)  -  u (•)  -  A  /  v?(u(o))do  in  P'(RN) 


This  la  easily  obtained  by  multiplying  (2)  by  test- functions  a  (t)9(x) 

n 

e  e  P(sF>  and  «n<t)  e  P(0,T)  converging  to  in  a  suitable  way. 

Note  that  tfc *  assumptions  on  <fi  together  with  ( 1 )  imply 


(5)  <?(u)  c  l’ut,!)  x  rN)  n  L**((t,t)  x  rN),  e  (0,T)  . 


Actually  the  relation  (4)  defines  u(t)  for  all  t  e  (0,T] .  Moreover  since 
*  IN 

u(t)  -  u(s)  and  J  *(u(o))do  e  L  (R  ),  for  all  0  <  s  4  t  4  T: 


(6)  /  u(t)  -  /  u(s) 

rn  rn 

t  t 

(7)  v(t)  -  v(s)  (A  /  <p(u(o))do)  »  -  /  i?(u<o))da  4  0  a.e.  on  rN  . 


In  (7)  we  use  a  uniqueness  result  (see  e.g.  [4]  lemma  A.5). 

The  relation  (6)  implies  the  relative  compactness  of  {u(t)  »  t  e  (0,T)} 
N  H 

in  c(M(R  ),  C  (R  )).  The  monotonicity  proves  the  uniqueness  of  the  limit  y 
of  u<t)  and  the  second  part  of  the  lemma  (see  e.g.  [17]  theorems  0.6,  3.8, 
3.9  about  potentials  of  measures).  Moreover 


v(t)  4  v(0)  **  /  -  Av(t)  4  /  -  Av(0)  . 

rn  rn 


Hence  /  u(t)  converges  to  J  -  Av(0)  -  f  dy  and  u(t)  converges 

*  „  I f  *" 

v  in  o(M(R  ),  CdR**)). 
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proof  of  theorem  1  for  N  >  3 

Let  h  e  ]0,T[  be  fixed.  If  u  and  u  are  solutions  of  (P),  by  (4)  we 

have 

(8)  VO  <  s  <  t  <  t+h  <  T, 

t 

(u(t)  -  u(t+h))  -  (u(s)  -  u(eth))  »  A  /  [^(u(a))  -  ^(u(o+h))]do  . 

s 


Letting  s  go  to  0  gives  in  O' (RN): 


t 

(9)  (u(t)  -  u(t+h))  -  (y  -  u(h))  ■  A  /  [^(u(o))  -  ip(u(o+h) )]  do  . 

0 

t  t 

Remark  that  s  f  ip(u(o))do  is  nondeoreasing  and  A  /  <p(u(o))do  is 

8  S 

bounded  in  L1 (R^) .  By  the  results  in  [4]  (lemma  A.5) ,  it  converges  in 

IN  r*' 

L.  (R  )  to  J  <i(u(o))do.  Let  us  denote 
roc  Q 

t 

g(t)  -  /  [<Mu(o))  -*(u(o+h)]do  +  v(h)  -  v(0)  . 

0 


Then  (9)  can  be  written  as 


u(t)  -  u(t+h)  *  Ag(t)  (■*— ^g( t)  ■  v(t+h)  -  v(t))  , 


where  v(t)  is  defined  in  lemma  2.  This  implies 


(10)  gfc(t)  «  a(t)  Ag(t)  , 


where 
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Tffr  Wiirinr  •  *  • 


a(t,x) 


g(u(t,x))  -  (u(t+h,x)) 


u(t,x)  -  u(t+h,x) 


if  u(t,x)  j*  u(t+h,x) 


if  u(t,x)  *  u(t+h,x) 


The  function  a  is  nonnegative  and  is  in  L  ((t,T)  x  R  )  for  any 

T  e  (0,T) •  Hence  g  is  solution  of  a  linear  equation;  moreover  if 

lim  u( t )  ■  lim  u(t),  g(0)  ■  v(h)  -  v(0)  is  nonpositive  by  lemma  2.  If 
t+  0  t+  0 

a(*)  were  regular  enough,  by  the  maximum  principle  applied  to  (10)  we  would 
obtain 


VO  <  t  <  t+h  <  T,  g(t)  -  v( t+h)  -  v(t)  <  0  . 


And  that  would  imply  v  <  v,  and,  by  a  symmetric  argument  v  ■  v  and  u  **  u. 
What  follows  is  to  justify  this  maximum  principle  for  the  equation  (10)  in  our 
particular  case.  The  method  consists  in  multiplying  (10)  by  the  solution  <(' 
of  the  dual  problem  +  A(ai|>)  ■  0,  ♦(T)  ■  0  e  J)+(RN),  0  <  T+h  <  T,  which 
formally  gives 


/  N  g<T)6  -  /  N  g(B)i>(B)  . 

R  R 


(V* (tF)  denotes  the  space  of  nonnegative  C  -functions  with  compact  support 

in  tF). 

Then,  we  show  that  the  right-hand  side  has  a  nonpositive  limit  when  s  +  0. 
The  first  step  is  to  "solve"  the  above  equation.  For  this  let  us 


approximate  a  by  a^  e  C  ( [0 ,T]  x  r  ) ,  nonnegative  and  satisfying: 


-6- 


N 

a  ,  |Va  |,  Aa  are  bounded  on  [0,T]  x  r  for  any  p 
P  P  P 

—  _  H 

vt  e  (0,7),  a  is  bounded  on  [t,T)  x  r  independently  of  p 
P 

_ 

a  converges  to  a  a.e.(t,x)  e  (0,T)  x  r  . 

P 

(For  instance,  one  can  mollify  a  and  multiply  by  a  C  -function  equal  to 


1  for  | x |  <  p  and  equal  to  0  for  |x|  >  p+1). 


Then,  e  >  0  being  fixed  we  consider  the  dual  problem 


(12)  |r  ♦  +  A(  (a  +c)«  )  -  0,  *  (T)  -  9  e  P+(RN)  , 

dt  p  P  P  P 


where  0  <  T+h  <  T.  For  simplicity  we  still  denote  T  by  T.  It  is  well- 

00 

known  that  this  problem  has  a  nonnegative  C  -solution  such  that 

i|>  (t)  e  L1  (RN)  n  L°(RN)  for  all  t  (see  [16]). 
p 

Now,  multiply  equation  (10)  by  i|>  and  integrate  to  obtain: 

P 

T  3i|> 

/  „  g(T)9  -  /  „  g(s)i|»  (8 )  ■  /  /  „  (g(o)  T“(o)  +  a(o)\|>  (a)Ag(o))do 

V  RN  p  s  RN  3t  P 

(13) 

T 

-  /  /  „  (a  -  a  -  e )t|»  (o)Ag(o)da  . 

8  rn  P  P 


In  order  to  pass  to  the  limit  in  p  for  s  e  (0,T),  let  us  make  some 
estimates  on  .  For  convenience  we  denote 


H  (t)  -  R  *  ♦  ft)  (-  -AH  (t)  “  ♦  ft))  . 
p  N  p  P  P 


Multiplying  (12)  by  Hp(t)  gives 


(14) 


I'*  1VHp(tH2  *  /  / 

R  R 


R 


(a  +e)^ 
P  P 
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This  proves  that  ♦  is  bounded  in  L2(0,T/RN)»  it  has  a  subsequence  weakly 

converging  to  <{f  .  Then  (a  +e)ij>  also  weakly  converges  to  (a+e)i|i  in 
e  p  p  t 

L2((t,T)  x  rN)  for  any  T  e  (0,T).  Hence  the  limit  \|>e  satisfies,  in 
an  integrated  form  of  (12),  namely 


VO  <  s  <  t  <  T  i}i  (t)  -  ib(s)  *  -A  /  (a+e)iji  . 

£  c  £ 

s 


Since  /  4  (t)  *  /  .  6,  ♦  (t)  is  bounded  in  L*(RN).  Hence  if>  is  in 

rN  p  rN  p  e 

L 1 ( RN )  and  by  (15) 


/  „  6  • 
R  R 


Moreover,  one  can  assume  that  iji  (t)  converges  to  i|>  (t)  in  o(M(RN) , 

p  E 


C^uF))  for  all  t  >  0. 


Now  (13)  becomes: 


A 

(17)  /  N  g(T)9  -  /  N  g(s)i|;e(s)  -  -e  /  /  R  ^(oJAgtOdo  . 


(Remember  that  Ag(t)  =  u(t)  -  u(t+h)  is  bounded  in  any  LP((s,T)  x  rN) 
and  Sp  converges  to  a  a.e.  on  (s,T)  x  R^  and  is  uniformly  bounded.) 

Now  we  let  e  ♦  0.  For  any  s  e  (0,T),  the  right-hand  side  of  (17)  is 
bounded  by 


elAgl  „  •  j  l*e(OI  1  N  , 

L  ( ( s , t )  X  r")  S  L  (R  ) 

00  IN 

which  converges  to  0  since  is  bounded  in  L  (0,T,I»  (R  ))  by  (16).  If 

H  (t)  -  E  *  (t),  integrating  (15)  gives 

£  No 


8 


- 


He(t)  -  H£{ a)  =  /  (a+e)ij»e 
s 


By  the  nonnegativity  of  t|/£  and  H£  this  implies 


0  <  H_(t)  <  H£(T)  =  *  6  . 


P  N  N 

Hence  H  is  bounded  in  L  ((0,T)  x  r  )  for  p  >  — —  and  one  can  find 
£  N-2 

convex  combinations  of  these  H£  converging  a.e.  and  strongly  to  H  in 

p  N  N 

L  ( (0,T)  x  r  )  for  some  p  e  (^^,»),  Since  ijj£(s)  is  uniformly  bounded 
in  L1^),  we  can  assume  that  the  same  combinations  of  <|/£( s)  converge 
a.e.  s  in  a(M(RN),  Cc(RN))  to  v(s)  e  M+(RN).  In  order  to  pass  to  the 
limit  in  (17),  we  need  a  convergence  in  o(M(R  ),  Cb(R  )).  This  comes  from 
the  fact  that 


Vs  e  ( 0  ,T) ,  /  dv(s)  =  /  9  -  lim  /  (s)  . 

N  #  N  _ _  £ 

R  R  £+0  R 


Indeed,  for  any  s  e  (0,T),  /  (a+e)i|i  is  bounded  in  L  ^  ( RN )  •  Hence,  by  (15) 

+  N  8 

there  exists  p(s)  e  M  (R  )  such  that 


0  -  v ( s )  =  -Ap(s)  in  V (RN) 


This  together  with  (16)  implies  (18). 

We  finally  obtain  the  existence  of  a  family  of  nonnegative  finite 
measures  {v(s),  s  e  (0,T)}  such  that 


(19)  Vs  e  (0 ,T)  /  N  g(T)0  =  /  N  g(s)dv(s) 

R  R 


-9 
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Li 


and 


(20) 


s  **■  H(s)  *  EN  *  v(s)  is  nondecreasing  on  (0,T) 


((20)  comes  from  the  monotonicity  of  H£). 


Now,  g(s)  *  v(s+h)  -  v(s)j  by  the  monotonicity  results  of  lemma  2,  (19) 


gives 


(21)  VO  <  s  <  3, 


/  N  g(T)0  <  /  (v(h)  -  v(s0))dv(s)  -  /  (u(h)  -  u(s  )>H(s>  . 

R  R  R 


(See  the  remarks  below  for  the  integration  by  part.)  But  j  dv(s)  is 

rn 

bounded  independently  of  s  and  H(s)  decreases  pointwise  when  s  decreases 
to  0.  Hence  v  *  -Ah(0+)  e  M+(RN)  and  by  (21) 


(22)  VO  <  sQ , 

/  N  g(T)0  <  /  N  (u(h)  -  u(Sp) )H(0+)  “  /  N  (v(h)  -  v(s0)) 


dv 


R 


Now  letting  s0  decrease  to  0  gives  by  monotonicity 


(23) 


/  „  9(T>e  *  /  M  <v(h)  -  v( 0 ) ) dv  . 

N  N 

R  R 


N 


But  v(0,x)  =  v(0,x)  >  v(h,x)  for  all  x  e  R  .  Hence 


V  e  t  V  (R  ),  /  g(T)0  <  0  . 


-10- 


A 


This  implies  the  relation  (11)  we  were  looking  for* 

Remark  1.  In  the  above  we  often  use  the  fact  that  given  p,v  e  Al+ ( RN ) : 


(24) 


/  (By  *  u>«  »>* 

R  R 


§ 


whatever  this  integral  is  finite  or  not.  In  (22),  H(0+)  is  the  decreasing 

limit  of  the  potentials  H(s).  It  is  generally  not  a  i.s.c.  potential  itself 
but  is  equal  a.e.  to  EN  *  v.  Since  u(h)  -  u(sQ)  is  a  "good"  function,  the 
integration  by  part  works.  It  would  not  for  h  m  0,  for  u(0)  is  only  a 
measure  (see  e.g.  [17]  for  more  details). 

Remark  2.  The  same  method  would  give  a  similar  uniqueness  result  for  the 
equation 


u  -  Avf(u)  in  V'  ((0,T)  *  0) 

<  Xu)  -  0  on  30"  (or  <p(u)  *  0  on  30") 

u(t)  +  d  in  o(M(RN>,  C.  vRN)) 

V  b 

with  0  a  regular  bounded  open  subset  of  R**,  by  using  the  "potentials" 
v(t)  solutions  of 


C  -  Av(  t)  “  u(t) 

\  v(t)  -  0  on  30  (or  v(t)  -  0  on  30)  . 

This  method  would  clearly  contain  the  cases  N  *  1*2. 

Remark  3.  There  is  no  potential  in  R1*  if  N  -  1,2.  Hence  we  have  to  do  the 


above  computations  in  an  "approximated"  way  using  the  solutions  of 


av  (t)  -  Av <  t>  “  u(  t )  , 
a  a 

and  latting  o  go  to  0 .  This  requires  more  a  priori  assumptions  on  the 

solution  u,  namely  1(u)  e  l/((0,T)  x  rN).  This  is  implied  by  (1),  (2)  in 

most  cases  of  interest,  like  the  porous  media  case  1(r)  -  rm  (see  theorem  4) 

.  In 

or  the  Stefan  problem  case  ^(r)  ■  (r-1)  (since  u  e  L  ((0,T)  x  r  )).  it 
was  also  proved  in  [14] ,  that  this  always  holds  in  dimension  1  with  very  weak 
assumptions  on  <p. 

Proof  of  theorem  1  for  N  °  1,2. 

Since  there  is  no  potential  in  R**  if  N  ~  1,2,  we  will  use  the 
solutions  of 

(25)  V  t  >  0  av  (t)  -  Av  (t)  -  u(t),  av  (t)  -  Av  (t)  ■  u(t)  . 

a  a  a  a 

For  a  >  0,  (a-A)  1  is  a  "good"  operator  even  when  N  -  1,2  and  in 
particular  (see  e.g.  [4]  lemma  1.1): 

u(t)  e  l1(RN)  n  l"(RN)  ■»  va(t)  t  tV>  n  C]t)(RN)  . 

We  will  denote  K  the  kernel  associated  with  (a-A)  \  i.e. 
a 

vQ(t)  -  Kq  *  u(t).  The  result  corresponding  to  lemma  2  is: 

Lemma  3.  If  u  >  0  is  a  solution  of  (P),  u(t)  converges  in 

o(M(RN)»  Cb(R^))  to  some  u  e  M+(RN)  when  t  ♦  0.  Moreover  when  t 

T 

decreases  to  0,  v  (t,x)  +  a  K  *  J  i(u(o))do  increases  to 
a  a 

t  t 

K  *  (y  +  a  /  v>(u(o))do(x)  for  all  x  e  rn  and  all  t  e  (0,T). 

°  0 


\ 
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T 

Mote  that  for  t  e  (0,T)  fixed,  K  *  /  ?(u(o))do  ie  oontinuoua  and 

“  t 

when  t  decreases  to  0,  it  increases  to  the  A. s.c.  function 

T  T 

K  *  /  ^(u(o))do  which  is  well-defined  since  /  v>(u(o))dc  e  L  (R  )  by 

°  0  0 

assumption  on  <fi. 


If  w  (t)  •  v  (t)  +  a  K  *  /  <Ku(0))do,  we  have  for  0  <  s  <  t: 

a  a  a  ' 


w  (t)  -  w  (s)  •  K  *(u(t)  -  u(s)  -  a  /  ^(u(0))do) 


t  t 

K  *((A-a)  f  l(u(0))do)  ■  -  /  ip<u(o))do  <  0  a.e. 
a 

s  s 


This  proves  the  second  part  of  the  lemma.  Then  we  finish  as  in  lemma  2  using 
that 


“  Vt!  *  ^„n  “  V01  ■  * 

R  R  R  R 


Now  to  replace  the  function  g  of  the  previous  proof,  we  introduce  for 
h  >  0  fixed: 


t  t 

g  (t)  ■  f  G(o ) do  +  v  (h)  -  v  (0)  -  o(a-A)”  /  G(o)do 

Cl  '  Cl  ot 

o  o 


where  we  denote  G{0)  «  ^(u(0))  -  ^(u(o+h)).  Then  we  verify 


(26)  (u-A)ga(t)  «  u(t+h)  -  u(t)  («-*  ga(t)  -  va(t+h)  -  v^Ct)) 


13 


(27) 


gat  +  «(a-A)ga(t)  -  -a(o-A)”1G(t)  <  a  *  *(u(t+h))  , 


where  a  is  defined  in  the  previous  proof* 

Now,  we  do  the  same  computations  as  for  N  >  3  with  the  operator  a-A 

instead  of  -A  and  we  obtain  -  like  in  (19),  (20)  -  the  existence  of 

{v  (s)  eM+(rN)>  s  e  (0 ,T)}  such  that,  for  6  eP+(R**)j 
01 


vs  e  (0  ,T) 

(28) 

T 

/  N  go(T)0  <  /  ga(s)dva(.)  +  /  /  a  *  «(u(o+h) )dwa(o)  , 

R  R  SR 


(29)  s  **•  H  (s)  -  K  *  v  (s)  is  nondecreasing  on  (0,T)  . 

<x  a  a 


By  lemma  3,  for  all  T  e  (0,T)  and  0  <  s  <  Sqi 


t 

g  (s)  <  [v  (h)  +  a  K  *  J  <p(u(o+h) )do) 
a  a  a  Q 

(30) 

T  T 

-  lv  (s„)  +  a  K  *  /  *(u(o))do]  +  a  K  *  /  G(a)do  . 

a  o  a  s.  a 

0  s 


Now  we  can  pass  to  the  limit  when  s  ♦  0  as  in  (22).  The  last  term  of  (30) 
is  easily  controlled  after  integration  by  part. 


i“  y  n  »«<•>*'.<»> 

S+0  R 

T  T 

<  /  Ntva(h)  +  a  X a  *  /  *(u(o+h))do  -  vo(sQ)  -  a  Kq  *  J  *(u(o) )do]dva 

R  a  a  0  sQ 


+  /  a  H  (0+)  /  G(0)do  , 

R  0 


14- 


with  v  -  li*  v  (a)  and  H  (0+)  ■  K  *  v  a. a.  .  Ha  lat  an  dacraasa  to 
®  s+0  8  °  ®  «  0 

0  above,  usa  tha  aonotonicity  established  in  leaaa  3  to  obtain 


lia  /  g  (s)dv  (s)  <  /  a  H  (0+)  /  *(u(o))<to  . 

a+0  *N  °  “  RN  °  0 


Finally,  cowing  back  to  (28)  and  remarking  that  H  (a)  <  H  (T)  -  K  *8  wa 

a  a  a 

have,  for  any  0  e  D+ (■?*), 


I j,  «.'*»• « I#  •  •  *.  * 


h  A  T 

/  *(u(o))do  +  J  I  a  8  K  *  <9  (u(o+h) )  . 

0  0  * 


Hanca 


T+h 

g  (T)  <  a  K  *  /  *(u(o) )do  . 

«  O  A 


Now  we  lat  h,  then  a  go  to  0* 
converges  to  0  in  P*  (■?)  whan 


For  any  (  t  L1 
a  ♦  0.  Indeed, 


(^*),  o 


r  •  a  X 

a  < 

*a  "  /  N  1 
H 


*  f 

iapliei 


tha  convergence  of  a  w_  to  son 

n  a 

n 


v  edV).  The  relation 


2 

aw  -  A  (a  w  )  ■  af  Implies  that  Av  -  0.  Hence  v  »  0.  Coming  back  to  tha 
a  a 

definition  of  g  (T),  we  obtain 
a 


T 

/  *(u(o))  -  l(u(o))  <  0  . 

0 


This  completes  the  proof* 


Remark*  Thanks  to  tha  laaaas  2  and  3,  tha  condition  (3)  in  theorem  1  can  be 
weakened  to  tha  requirement 


15* 


(3)' 


11*7/  (u(t)  -  u<t))dt  -  0  in  PM**)  . 

h+0  h  0 

This  any  be  uatfnl  in  view  of  the  existing  literature  where  the  solutions  ere 
often  defined  as  functions  u  satisfying  (1)  and 

T 

/  /  u  u*  +  *(u)A*  +  /  M  u(0 )♦(<>)  -  0 

0  K*  t  IT 

«  U  A 

for  any  f  e  c  (tO,T[  *  R  )•  Clearly  two  solutions  u  and  u  of  the  above 
satisfy  (2)  and  (3)'. 
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initial  data, 
exists  a  unique  non- 
all  t  >  0  such  that 

(31)  ufc  -  Aum  in  V  * ( (0,“)  x  RN) 

(32)  u(t)  ♦  p  in  o(M(RN),  C.  (RN))  when  t  +  0  . 

D 

A  +  N 

If  u  is  another  such  solution  with  initial  data  p  e  M  (R  ) 


SECTION  II.  Existence  results  and  dependence  on  the 

+  N 

THEOREM  4.  Let  m  >  1  and  p  e  M  ( R  ) .  Then  there 
negative  u  e  C((0,")j  L1 (RN ) )  n  L  ((t,«)  x  rN)  for 


(33) 


Vt  e  (0,»)  /  N  |u(t)  -  u(t)|  <  /  N  Iw-pI  . 


Moreover,  if  p  e  M+(RN)  converges  to  p  in  o(M(RN),  Ch(l^)),  the 
n 

associated  solutions  u^t)  converge  to  u(t)  in  L ^  ( RN )  for  all  t  >  0. 

Rantark  4.  If  p  is  the  Dirac  mass  6  at  the  origin,  the  solution  of  (31), 
(32)  has  been  explicitly  determined  (see  Barenblatt  [2]).  It  is  given  by 


u(t,x) 


k(m-1 )x 


2mN  t 


k/N 


where  k-1 
way  that 


■1  + 


°(t’ 


2 

N 


and 

1. 


a 


is  a  constant  depending  on  m  and  N 


in  such  a 


Remark  5.  The  proof  of  the  above  result  contains  several  ingredients.  First 
the  existence  of  a  solution  to  (31)  when  the  initial  data  is  regular.  This 
was  proved  in  [20].  It  can  also  be  obtained  as  a  consequence  of  the  abstract 
theory  of  evolution  equations  governed  by  accretive  operators  which  carries 
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over  to  the  more  general  equation  (2)  (see  [5]).  Then,  approximating  w  by 

"regular"  functions  u  ,  one  has  to  prove  that  the  solution  of  (31)  with 

n 

initial  data  w  converges  to  the  solution  of  (31)  -  (32).  This  needs  some 
n 

compactness  arguments  that  can  be  obtained  through  two  different  ways.  One 
can  use  that 


I  Aum| 


C  t 


as  proved  in  [1].  Actually  this  method  would  apply  to  (2)  for  the  general 
class  of  functions  <fi  described  in  (8).  Here  we  will  use,  in  conjunction 

m 

with  som  direct  estimates  in  (31),  the  L  -regularising  effect  which  says 
that  the  solution  of  (31)  belongs  to  L**((t,")  x  r*1)  for  any  T  >  0.  The 
latter  property  -  which  is  needed  to  apply  our  uniqueness  result  -  is  also 
true  for  equation  (2)  with  very  weak  assumptions  on  <P»  To  illustrate  the 
generality  of  this  method  let  us  establish,  at  least  for  N  >  3,  a  more 
general  existence  result. 

Let  us  consider  i  [0,**)  ♦  (0,**)  locally  Lipschitz,  nondecreasing, 

*(0)  >0.  Assume  N  >  3  and  for  instance  (see  [3])s 


(34) 


3a  > 


N-2 

N 


such  that  (ip  ( r) ) 


1/a 


is  convex  for  r  large 


Then  we  have 


PROPOSITION  5.  Por  all  it  eM+{RN),  there  exists  a  unique  non¬ 
negative  u  e  l"((0,»)»  L^R*1))  n  x  rN)  for  all  t  >  0,  solution  of 

ufc  ■  Aip(u)  in  P'  ( (0,*)  x  rN) 

u(t)  ♦  y  in  o(M(rN),  C^(RN))  when  t  +  0  . 

Moreover  the  estimate  (33)  holds* 

00  M 

Remark  6.  The  assumption  (34)  insures  that  uel((t,**)xr)  for  t  >  0. 

A  slightly  different  assumption  can  be  found  in  [21]* 

Proof  of  Proposition  5. 

In  n 

Let  wn  e  L  (R  )  nonnegative  and  converging  to  y  in  o(M(R  ), 

cb(R«>).  By  the  existence  results  in  [5]  and  the  L**-regularising  effects 

established  in  [3],  there  exists  u  e  C([0,»)i  L^R1*))  n  l"((t,«)  x  r**)  for 

n 

any  t  >  0  such  that 

(35)  u  **  A*(u  )  in  P’((0,<»)  x  rn) 

nt  n 

y 

,  K,  KQ,  Y  >  0  , 

where  K,  KQ ,  Y  depend  only  on  ly  I  ^  o,  t p  and  N.  Remark  that  is 

a  »  L  ^ 

uniformly  bounded  in  C((0,»)»  l  (R  ))  since  /  „  u  (t)  ■  J  „  y  by  (35). 

_N  n  _N  n 

R  R 

Let  us  make  some  formal  estimates  now.  Multiplying  (35)  by  <P(un)t  yields: 


(36) 


lu  <t)l  <  K  + 
"  L 


(37) 


''•  fc  ••  •-*■•  . 


!P9P«^P£pipp!rr 


*<£  .  -  .... 


’4 It  I  »  l,',«“n"2  ’  *(jl  V  "'“n’t  ‘  0 


R 


Multiplying  (35)  by  <*>(un)  and  integrating  give  for  0  <  t  <  t 


(38) 


f  J  |V*(un)|2  -  /  N  4-(un(T))  -  i(i(un(t))  , 
T  RN  " 


where  <l»(r)  ■  f  ^(o)dtj.  Since  u  (t)  is  bounded  in  I>  independently  of  n 
0  n 


(see  (36)),  this  implies  that  V<tfu  )  is  uniformly  bounded  in 


L2  ( (0 ,«•)  x  RN).  By  (37),  7^(u  )  is  even  uniformly  bounded  in 

loc  n 


lT{x,mt  L2(RN))  for  any  t  >  0.  Hence,  integrating  (37)  proves  that 

OS 

f  f  u  <p(u  )  is  bounded  for  any  T  >  0.  Since  V  is  locally 
J  ‘  _N  nt  n  t 


T  R 


Iiipschitz ,  we  obtain  that  /  /  N  [V»(un>tJ2  is  bounded.  Finally  we  deduce 

T  R 


that 


uj(u  )  is  bounded  in  wl'^tO,0")  x  rN) 
n  roc 


(39) 


the  bound  depends  only  on  sup  1 Un ■ ^ 


The  formal  computation  (38)  can  be  justified  like  in  [5]  -  Prop.  10.  For  the 
other  ones,  we  use  <p(u(t+h))  -  <p(u(t))  and  let  h  go  to  0. 

By  (39),  there  exists  a  subsequence  (still  denoted  ¥>( un) )  converging  in 


l?  ((o,<»)  x  rN)  and  a.e.  to  some  w.  On  the  other  hand,  a  subsequence  of 
toe 


un  converges  weakly  in  IT(K)  for  any  compact  K  c  (0,»)  x  r  and  the 


limit  u  satisfies  w(t,x)  -  <?(u(t,x))  a.e.  since  *  is  a  maximal  monotone 
operator  in  L2(K)  (see  [7]  prop.  2.5).  Clearly 


u  e  L*(0,T,  L1  (*?*))  n  l*((t,«)  x  rN)  for  all  T  >  0  and  satisfies 
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ut  *  A^(u)  in  V'  ((0,®»)  x  RN)  . 

It  remains  to  show  that  0  *  lim  ess  u(t)  (in  o(M(RN),  Cb(RN))  which  exists 

t+0  f  A  f 

by  lemma  1 ,  is  equal  to  y  ( note  that  J  dy  <  J  dy ) . 

For  this,  let  us  assume  we  have  chosen  y  «*  y  *  p  with 

n  n 

M  M  00 

p  (x)  »  X  n  p(n|x|)  fx  e  R  where  pec  ([0,«»)]  is  supported  in  [0,1]  and 
n  N 

X„  is  a  constant  such  that  /  „  p  *  1.  with  this  choice  (see  [17]) 

N  _N  n 

R 

n  *  vn  *  *  Wn  *8  nondeoreasin9  (and  v°  increases  pointwise  to 

v°  *  E  *  y ) •  Hence 


n  ♦  v  (t)  ■  E  *  u  (t)  is  nondecreasing  . 
n  N  n 

Indeed  v„(t)  is  a  solution  of  v  .  +  <p(-&v  )  *  0,  v  (0)  -  v°  and  one  can 
n  nt  n  n  n 

use  the  maximum  principle  for  this  equation  (see  e.g.  [9]).  Since 

/  „  -  Av  (t)  ■  /  „  u  (t)  is  bounded,  v_(t)  increases  to  a  potential  v(t) 
J  n  rN  n  n 

such  that  -Av(t)  is  the  limit  (in  o(M(RN),  C_(RN))  of  -Av  (t)  (see  e.g. 

c  n 

[17]).  Necessarily  v(t)  «  Ej,  *  u(t)  (at  least  a.e.  t).  Now  ,by  lemma  2,  if 

A  ©  A 

v  “  E *  y,  we  have 
N 


V  >  V(t)  >  V  (t) 
n 


Vn,  a«ei 


-  v°  > 


-0  N  0 
V  >  V 


on  the  other  hand 


v°  >  v°  >  v  (t)  Vn,  a.e.  t 
n  n 

_  0  v  0*0. 

*  V  >  v(  t)  V  >  V 


21 


Hence  vQ  ■  vQ  and  y  =  y. 

To  complete  the  proof,  let  us  prove  (33).  By  accretivity  in  L^(R**) 
(see  [5]),  for  any  n,  t: 


/  N  |un(t)  "  \(t)l  <  ^  N  K  -  Kl 

R  R 


/  N  |Pn  *  (w-p ) I  <  /  N  Iu-pI  • 

R  R 


We  apply  a  Fatou-type  lemma  to  finish. 


Proof  of  theorem  4. 

For  N  >  3,  the  existence  of  u  is  a  consequence  of  the  proposition 
3.5.  Using  the  particular  structure  of  <P(u)  ■  uF,  we  add  an  argument  to  the 
previous  proof  in  order  to  absorb  the  cases  N  •»  1,2  and  to  prove  the 
continuity  results  for  all  N. 

Let  y  and  u  defined  as  in  this  proof.  The  estimates  we  established 
n  11 

are  valid  for  any  N.  Hence  (39)  holds  and  a  subsequence  of  converges 

a.e.  to  u  e  l”((t,~)  x  RN)  n  L°°(0,T,  L1  (R** ) )  solution  of  (31).  Moreover,  in 
this  particular  case,  we  have 

T 

lim  /  /  um(a)du  =  0  uniformly  in  n  . 

4  4  N  n 

T+0  OR 

Indeed,  by  the  L°°-estimate  (see  [3],  [21]) 


lu  (t)«  <  ~  In  I6, 

n  ~  n  1 


with  o 


m~1+i 


,  S 


2k 

N 


Hence 


22- 


* 


/•T  r  m,  .  .  „  „1  +  5(mr-1)  fT  dt 

L  Un(0,d°  <  C,VT1  Jn 

OR  L  0  t 


where  (m-1 )o  <  1 


In  particular,  for  all  t,  /  um(o)do  e  L(RN)  and  is  the  limit  in 

0 

l]  (RN)  of  /  um(a)do.  By  lemmas  2  and  3,  u(t)  converges  in 
roc  0  n 

o(M(RN),  C.  (RN))  to  some  jj  when  t  +  0.  Now,  we  can  pass  to  the  limit  in 


c  N 

u  (t)  -  p  =  A  /  um ( a ) do  in  P'(R  ) 


and  obtain  that  m  =  u. 

Remark  that  /  =  “(*>•  Hence  uR(t)  cc 

rN  n  R  R  R 

verges  to  u(t)  in  L  ^  ( R*1 )  for  a.e.  t  and  even  Vt  >  0  by  the 
contraction  property.  The  uniqueness  proves  the  convergence  of  the  whole 
sequence . 

N 

By  (39)  for  any  open  subset  H  relatively  compact  in  (0,«°)  x  R 


*u“l  12  <  C(^  NU)  * 

W  '  <fl)  R 


And  (41)  gives 


/  J  um(0)d O  <  C  •  T 

0  rn 


1-o (m-1 ) 


[/rN  .] 


1+5 (m-1 ) 


For  the  uniqueness  when  N  =  1,2,  given  u,  u  solutions  of  (31),  (32) 
we  apply  theorem  1  to  u(*  +  t)  (as  well  as  u(*  +  t))  for  any  x  >  0  to 
prove  that  they  coincide  with  the  solutions  in  the  sense  of  semigroups.  Hence 
(43)  holds  for  u,  u  and  we  can  apply  the  theorem  1  to  u  and  u. 


-23- 


M  .e 

Now,  let  UR  converge  to  u  in  o(M(R  ),  C^(RW))  and  u^,  u  the 
corresponding  solutions.  The  same  arguments  as  above  using  the  estimates 
(42),  (43)  plus  the  uniqueness  result  prove  that  uR(t)  converges  to  u(t) 
in  L1  (R**)  . 
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